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Abstract 

We study the CP violating asymmetry, the forward-backward asymmetry of the lepton pair and 
the CP asymmetry in the forward backward asymmetry for the exclusive decay B K*t^t^ 
in the general two Higgs doublet model including the neutral Higgs boson effects. We analyse 
the dependencies of these quantities on the model III parameters, We found that the physical 
parameters studied above are at the order of the magnitude 1% and neutral Higgs boson effects 
are detectable for large values of the coupling 
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1 Introduction 



Rare B decays are induced by flavor changing neutral currents (FCNC) at loop level in the Standard 
model (SM) and they are rich phenomenologically. They open a window for the determination of free 
parameters in the SM and also the investigation of the physics beyond, such as, two Higgs Doublet 
model (2HDM), Minimal Supersymmetric extension of the SM (MSSM) (jH], etc. The experimental 
work for rare B decays continue at SLAC (BaBar), KEK (BELLE), B-Factories, DESY (HERA-B) 
and this stimulates the theoretical effort on them. 

The exclusive B K*l^l^ process is an important candidate among rare B decays. Since its 
branching ratio (Br) predicted in the SM is large, there is a strong hope that this physical quantity is 
measured in the near future. The inclusive decay which induces B K*l^l^ process is 6 ^ sl^l~ 
transition and it is extensively studied in the literature, in the framework of the SM, 2HDM and 
MSSM. In P]-[[T5|], b — > sl^l^ process is studied for light lepton pairs, namely I = jj. In this case, 
the neutral Higgs boson (NHB) effects can be neglected since those contributions are proportional to 
the light lepton masses or corresponding Yukawa couplings. However for / = r, the NHB effects 
give sizable contributions. In [ [T^ , [TT]] B — > X^r+r^ process was studied in the model I, II versions 
the 2HDM and it was shown that NHB effects are important for large values of tan/3. Currently the 



inclusive b — > sl^l~ decay was studied in the model III version of 2HDM [[TS]] and it was observed 
that NHB effects can give considerable contribution if the Yukawa interaction between r lepton and 
neutral Higgs bosons is large. 

The theoretical analysis of exclusive decays is difficult due to the hadronic form factors which 
contain uncertainities. However, their experimental investigation is easier compared to the one for the 
inclusive decays. The calculation of physical observables in the hadronic level needs non-perturbative 
methods. In the literature there are different studies based on different approches, such as relativistic 
quark model by lightfront formalism [[T5|], chiral theory [|T9|], three point QCD sum rules method [pO|], 
effective heavy quark theory [ pi] ] and light cone QCD sum rules [ p2| , p3| ]. 

With the measured upper limit 5.2 x 10"^ (4.0 x 10~^) for the Br of the decay 5+ K~^jj,~^fi~ 
(5° — > K*^fi'^ji^) [^, the process B — > K*l^l^ have reached great interest. There are various 
studies on these decays in the SM, SM with fourth generation, multi Higgs doublet models, MSSM 
and in a model independent way, in the literature @]-[jl5l] and [lT9|]-[]36[]. 

The CP violating effect is an important physical quantity to ensure the information about the free 
parameters of the model used. Since the CP violation for B K*l^l^ decay almost vanishes in 
the SM due to the unitarity of Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and smallness 
of the term VubV*^, we have a chance to investigate the physics beyond the SM by searching the 
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CP violating effects. In the model III version of the 2HDM, there is a new source for CP violation, 
namely complex Yukawa couplings. In [ p6| , CP violating effects due to new phases in the model 
III and three Higgs doublet model, 3HDM(02), were studied and it was observed that a considerable 
CP asymmetry was obtained. 

In this work, we study the exclusive B K*t^t^ decay in the general 2HDM , so-called model 
III, by including NHB effects. We use the quark level effective Hamiltonian which is calculated 
in [ JTSl ] and investigate the Acp and the forward-backward asymmetry (Aps) of the lepton pair for 
the process underconsideration. Further, we calculate the CP asymmetry in Aps, {AcpiApp)) and 
observe that it can be measured in the forthcoming experiments. 

The paper is organized as follows: In Section 2, we present the leading order (LO) QCD corrected 
effective Hamiltonian and the corresponding matrix element for the inclusive b — > sr+r^ decay, 
including NHB effects. Further, we give the matrix element for the exclusive B K*t^t^ decay 
and the explicit expressions for ApB, Acp and Acp{Apb)- Section 3 is devoted to the analysis of the 
dependencies of App, Acp and Acp{Apb) on the CP parameter sinO, Yukawa coupling ^^^^ and 
the mass ratio and to the discussion of our results. In Appendices, we give the explicit forms of 
the operators appearing in the effective Hamiltonian, the corresponding Wilson coefficients and the 
form factors existing in the hadronic matrix elements. 

2 The exclusive B K*t^t~ decay in the model III including 
NHB effects. 

In the model III, the flavour changing neutral currents in the tree level are permitted and various 
new parameters, such as Yukawa couplings, masses of new Higgs bosons, exist. Yukawa couplings 
describe the interaction of fermions with gauge bosons. Our starting point is the inclusive b sr^r^ 
process which induces the exclusive B K*t^t^ decay and the corresponding Yukawa interaction 
is given by 

+ VkMiEiR + ^j^JkL(p2EiR + h.c. , (1) 

where i,j (k, I) are family indices of quarks (leptons), L and R denote chiral projections L{R) = 
1/2(1 =1= 75), (f)m for m = 1,2, are the two scalar doublets, Qn (hp) are quark (lepton) doublets, 
UjR, Djp (Eip) are the corresponding quark (lepton) singlets, rjfj'^ and ^fj'^ are the matrices of the 



2 



Yukawa couplings which have complex entries in general. Here (pi and 02 are chosen as 



1 

71 



\ ^( V2x' 



J_( V2H+ 
v/2 Hi + tH2 



(2) 



with the vacuum expectation values, 



^(0 



< 01 >= -= ( 1 ; < 02 >= . (3) 



V 



With this choice, the SM particles can be collected in the first doublet and the new particles in the 
second one. Further, we take Hi, H2 as the mass eigenstates h^, respectively. Note that, at tree 
level, there is no mixing among CP even neutral Higgs bosons, namely the SM one, H'^, and beyond. 

The part which produce FCNC at tree level is 



Cy,FC = ^fj^Q^L<p2UjR + egQ.L02/^j7J + iklhL(p2Em + h.C. . (4) 

In eq.(Q) the couplings ^^'^ for the PC charged interactions are 



^ci ~ ^neutral VcKM , 

ich = VcKM ineutral , (5) 



where Cneutrai is defined by the expression 

^U{D) _ /T.^7(D)^_l^{/,(D)T.C/(^)) .^N 
?iV — {^R{L) ) S ^LiR) ■ W 

and ^neutral dcnotcd as Here the charged couplings are the linear combinations of neutral 

couplings multiplied by Vckm matrix elements (see p7[ ] for details). 

At this stage we would like to present the calculation of the matrix element for the inclusive 
b ST^T^ decay, briefly. The procedure is the following: 

• The calculation of the full theory including the NHB effects which comes from the interactions 
of neutral Higgs bosons /i° and with r lepton. 

• Overcoming the logarithmic divergences by using the on-shell renormalization scheme. Here 
the renormalized vertex function is taken as 

yRen r 2\ ^ -pO / 2\ , yC /yx 
neutr \f J ^ neutr yf J ' ^ neutr i V / 

with the renormalization condition 



and the counter terms are obtained. Here the phrase neutr denotes the neutral Higgs bosons, 
H^, and and p is the momentum transfer. Note that the self energy diagrams do not 
contribute in this scheme. 

• Integrating out the heavy degrees of freedom, namely t quark, W^^ iJ^, hP, and bosons 
in the present case and obtaining the effective theory. 

• Performing the QCD corrections through matching the full theory with the effective low energy 
one at the high scale /i = mw and evaluating the Wilson coefficients from my/ down to the 
lower scale ^ ~ 0{mi,). 

• Obtaining the effective Hamiltonian relevant for the process b sr^r^ which is given by 

where Oi are current-current (i = 1,2), penguin (i = 3, ...,6), magnetic penguin (i = 7,8) 
and semileptonic {i = 9, 10) operators. Here, Cj(/i) are Wilson coefficients normalized at the 
scale /i and given in Appendix B. The additional operators Qi{i = 1, .., 10) are due to the NHB 
exchange diagrams and Cq.{iA are their Wilson coefficients (see Appendices A and B) . 

Therefore the QCD corrected amplitude for the inclusive h — > sr+r~ decay in the model III reads 

as, 

M = ^^V,,r;ia^^(s7^P,.6)f7^r + C'io(5-7;.i^L&)r7^75r 



7r 



- 2C',ff'^{sia^,q,PRh)f^^T + CQ,{sPRh)fr + CQ,{sPRh)f^,T^ . (10) 

The matrix element for B — > K*l^l^ decay can be obtained by inserting the inclusive level 
effective Hamiltonian in eq. (^ between inital, B, and final, K* , hadronic states. The neces- 
sary matrix elements in this calculation are {K* |s7^t(l ± 75)&| B), {K* Isiaf^^q^^l + 75)6! B) and 
{K* |s(l ± 75)61 B). They are calculated by using some non-pertubative methods like QCD sum 



rules, light-cone QCD sum rules, etc., and using the parametrization of the form factors as in [20], the 



matrix element of the B K*t^t decay is obtained as p2|]: 

Ga 
2V2n 



M = -^^V,bV*lfYr \2Ae,,p„e*^p'j,,q'' + iB^e* - iB^{e*q){pB + Pk*)^. - iB;{e*q)q^ 



2Ce^.p.e*X.g" + iDie^ - iD2{e*q){pB + Pk*)^. - iD^{e*q)q^ 



+ ifTF{t*q)+if^^TG{e*q) \ , (11) 



where e*'^ is the polarization vector of K* meson, pb and pk* are four momentum vectors of B and 
K* mesons, q — Pb —px*- A, C, F and G, Bi and Di, i — 1, 2, 3 are functions of Wilson coefficients 
and form factors of the relevant process. Their explicit forms are given in Appendix C. 

Now we are ready to calculate the forward-backward asymmetry of lepton pair, CP-violating 
asymmetry and CP violating asymmetry in forward-backward asymmetry for the given process. 

The forward-backward asymmetry ApB of the lepton pair is a measurable physical quantity which 
provides important clues to test the theoretical models used. Using the definition of differential ApB 

with z — cos 9, where 9 is the angle between the momentum of B-meson and that of r" in the center 
of mass frame of the dileptons t+t~, we get 

JdsE{s) 



S dsD{s) 
Here, 
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E{s) = QmBXv^i r T-AAmBmrRe{C',^^*CQMm,{^-l)A2{q^) 

yrribr [r - 1) \ ^ 

+ m6(v^ + 1) Ai {q^) -2mBsTs ((r - 1) (3 r - s + 1) T2 {q^) + (r' + (s - if 

- 2r(s + l))r3(g2 

2 (\ ("^B Re{Cl"* J (m,(^/f - 1) A^ {q') + m,{^ + 1) Ar {q') 

- 2mBsn {{l + V^fir + s- 1) A, (g^) + XA2 {q^' 
+ 8 Cio (| - 2 ms m;, ( - 1) Re {Cff) {q') V (g^) 

+ A^ {q') (2 m, mB{V^ + 1) Re {C',^^) T, {q') + m| s ReiC^") V{q'))^ | (14) 

D{s) = VXvl^-^^ml\Cr\H2ml + mls)(^--^ 

+ 8 T?(q^) A + + I n{q^) (12 (r - 1) V - (4 r - .) A)) 

+ (1 + J^)2^, ^^l^9^^l'(2^r + < s) (2 A,{q') A,{q') (1 + VF)^ (r + s - 1) A 

+ Al{q^) (1 + V^)^ (12 r s + A) + A (8 r s V\q^) + Al{q^) \)\ 



+ 2 Ci'o ms l^-i- (^2 A^iq^) A2 (g') (2 m^(r - 2 s - 1) + s (r + s - 1))A) 

—hAmBml{A2{q') (-1 + V^) + A^{q^) (1 + V^)) Tglg^) a) 

777,6 r 

+ , \o — ^2(9') f^il, s A + 2 rr4(6 s (1 + r) - 3 + A 



(1 + ^JrYr s 

+ ^ ((1 + V^' Al{q^){ml^ s(12 r s + A) + m^(-48 r s + 2 A)))^ 

+ |C7c?, r K. « - 4 m^) (A,{q') m, ( - 1) + ^(g^) (v^ + 1) - 2 T,{q') rriBs)' 

+ |c^Q.r (V^ - 1) + A(a (V^ + 1) - 2 r3(g2) V^. 

7775 r ^ " / 

+ n ^ V^2 ^e(C,^^^* Cg^^Oie (2 ml + s) (8(1 + ^) Uq^) V{q') A 

- — r^^U(g^) (A(r - 1)(1 + 3r - T^lg^) + X^q^)) + ^^(g^) (1 + 

( (r - 1) T2(g^) (12(r - 1) r - A) + (r + s - 1) Uq') A))) 

Cio i?e(CQj m,(m6 ((-1 + v^) A^iq^) + (1 + A,{q^)) 



ml r 



- 2mBn{q'))(A2{q'){l - V^) - A,{q') (1 + V^) + 2m|^ ^^^^3(52)) a| (15) 

where A = 1 + + — 2r — 2s — 2rs, r = ^^^^ and s = -4-. 

The NHB effects bring new contribution to and we will study those contributions in the 
Discusssion part. 

The complex Yukawa couplings are the possible source of CP violation in the model ni. In our 
calculations we neglect all the Yukawa couplings, except ^6 ^n,tt choose ^^^^ com- 

plex, — I ^^ 66 1 (s^^ Discussion part). Therefore the CP violation comes from the Wilson 
coefficients 6*7^^, Cg^ and Cq^. Using the definition of Acp 

_ TjB ^ K*T^T-) - T{B ^ K*T+T-) 

r{B K*r+r-) + r{B K*t+t-) ' ^ ^ 



we get 



where 



1 - ^dsm 



A A^{q^) ((r - 1)(3 r - s + 1) T^{q^) + m{q^ 



r{^/r - 1) 

+ (Vf + If A,{q') ((r - 1) T2(g2)(12 (r - 1) r - A) + (r + s - 1) A T3(g2))^ | 

and 



(18) 



Ais) = Dis) + Dcpis) . (19) 
Here Dcp{s) is the CP conjugate of D[s) which is defined as 

Dcp{s) = D{s){i^^,,-^i^%) . (20) 

The CP violating asymmetry in Air^ is also a measurable physical quantity and it can give strong 
clues for the physics beyond the SM. This quantity is defined as 

Acp(A,s) = f"'//' ■ (21) 

/ipB + ^FB 

where Afb is the CP conjugate of Afb and it is given as 

ApB = AFBi^N,bb ~^ ^N*bb) ■ (22) 

Note that during the calculations of Acp, ApB and Acp{Afb) we take into account only the 
second resonance for the LD effects coming from the reaction h s^i st^t^ , where i = 1, .., 6 
and divide the integration region for s into two parts : < s < and (™'^2^°-°^) < s < 1, 

where m^^ = 3.686 GeV is the mass of the second resonance (see Appendix B for LD contributions). 

3 Discussion 

In the general 2HDM model, the number of free parameters, namely the masses of charged and neutral 
Higgs bosons and complex Yukawa couplings increases compared to the ones in the SM and 

model I (II) version of 2HDM. The arbitrariness of the numerical values of these parameters can be 
removed by using the restrictions coming from the experimental measurements. 

Since the neutral Higgs bosons, and A^, can give a large contribution to the coefficient Cj^^ 
(see the Appendix of [RBp for details) which is in contradiction with the CLEO data p^]. 



Br{B X,7) = (3.15 ± 0.35 ± 0.32) 10"^ , (23) 
we take ^^^^ ~ and ~ 0, where the indices i,j denote d and s quarks . Further we use the 



constraints p7| 1, coming from the AF = 2 mixing (here F = K, Bd, D) decays, the p parameter [^011 , 



and the measurement by CLEO Collaboration eq. (|23|), we get the condition for ^Tvtc. ^Ntc « (.mt 



and take into account only the Yukawa couplings of quarks C,N,tt ^i^d C,N,bb- We keep the Yukawa 
coupling ^^^^ free and increase this parameter to enhance the effects of neutral Higgs bosons. 

hi this section, we study the CP parameter sinO, the Yukawa coupling and the mass ratio 



1^0 



dependencies of the ApB, Acp and Acp{Apb) of the exclusive decay B K*t^t , restricting 



\Cf^\ in the region 0.257 < \Cj^^\ < 0.439 due to the CLEO measurement, eq.(|2|) (see ^ for 
details). Our numerical calculations based on this restriction and throughout these calculations, we 
use the redefinition 



U,D 



lAG 



V2 



we take I 



< 1, the scale /i = m^, include the LD effects and use the input values given in Table 



Parameter 


Value 


TTLr 


1.78 (GeV) 




1.4 (GeV) 


rUb 


4.8 (GeV) 


P^ 


40 rrib 


"em 


129 


Xt 


0.04 




175 (GeV) 


mw 


80.26 (GeV) 


mz 


91.19 (GeV) 


rriHO 


150 (GeV) 


rriho 


70 (GeV) 


mH± 


400 (GeV) 


^QCD 


0.225 (GeV) 


as{mz) 


0.117 


sinOw 


0.2325 



Table 1: The values of the input parameters used in the numerical calculations. 
In Fig. [D we present sinO dependence of App without NHB effects, for m^o = 80 GeV. Here 



ApB lies in the region bounded by solid (dashed) lines for C^-^-^ > (Cj^^ < 0). The solid straight 
line shows the SM contribution. In the model III without NHB effects, \App\ is smaller compared 
to the one in the SM (0.195), for Cj-^'^ > 0, however it is possible to enhance it at the order of the 
magnitude 2% with increasing sinO. For Cj-^^ < 0, App is not sensitive to sin6 and the restriction 
region is narrow. For this case \App\ can have slightly greater values compared to the SM one. 
Addition of NHB effects (see Fig. ^ reduces \Apb\ for C^-^-^ > almost 30% compared to the one 
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without NHB effects. For C^-^-^ < 0, the restriction region becomes narrow and ApB reaches the SM 
prediction for small sin9. 

Fig. D represent Cw.rr dependence of ApB for sinO = 0.5 and m^o = SOGeV. \Apb\ vanishes 
with increasing ^^^^ for C^-^-^ > 0. For Cj-^'^ < 0, \ Apb\ does not vanish in the given region of ^^^^^ 
and it stands less than the SM result. 

Fig. is devoted to the ratio dependence of ApB for sinO = 0.5 and ^^tt — 10 m^-. 
Increasing values of the ratio causes to increase \Afb \ for both C^-^^ > and C^^-^ < 0. If the masses 
of and A^ are far from the degeneracy, \Afb \ becomes small especially for Cy"^-^ > 0. 

Figs HQ represent Acp of the process B K*t^t^ . In Fig. |] we present sinO dependence of 
Acp without NHB effects, for wlao = 80 GeV. Here Acp lies in the region bounded by solid (dashed) 
lines for C^^^ > (Cj^^ < 0). For Cj^^ > 0, Aqp is at the order of the magnitude of 1% for the 
intermediate values of sinO and its sign does not change in the restriction region. However Acp can 
have both signs, even vanish for Cj-^-^ < 0. With the addition of NHB effects (see Fig. ||) Acp for 
Cj^^ > decreases to almost one half of the value we get without NHB effects. For Cj^^ < there 



is still a decrease in Acp.This behavior can be seen from the expression eq. (jTSJ) since the numerator 
of the Acp ratio is free from NHB effects and their additional contributions enter into the expression 
in the denominator part. Further the restriction regions becomes narrow. 

Fig. represent ^^^^ dependence of Acp for sinO = 0.5 and m^o = 80 GeV. Acp is sensitive 
to the parameter .^^^^ and it decreases with increasing ^^^.^ for C^^-^ > 0. However, for Cj-^^ < 0, 
the dependence of Acp to ^^^^ is weak. 



The ratio dependence of Acp for sinO = 0.5 and C,^^-^ = 10 m^- is presented in Fig. |[ As 



seen from the figure the sensitivity Acp to the ratio is small, especially for Cj < 0. 

Finally, we present the CP violating asymmetry in in a series of figures (Figs. P[|n|). Fig. ^ 
represent sin9 dependence of Acp{Afb) with NHB effects, for m^o = 80 GeV and ^^tt = 10 "^t- 
Acp{Afb) is at the order of the magnitude of 1% for the intermediate values of sin9 for G^^-^ > 0. Its 
sign does not change in the restriction region similar to the Acp of the process under consideration. 
However Acp{Afb) can have both signs, even vanish for Gj^^ < 0. Acp{Afb) is sensitive to the 
parameter .C^,rr especially for the large values of ^^^^t ^^'^ ^i^^ > (see Fig. |T0|). It can reach 10% 







for ^^^^ = 50 GeV . In the case Gj^^ < 0, Acp{Afb) is not sensitive to Cwrr almost vanishes. 

Fig. [TTJis devoted to the ratio dependence of Acp{Afb) for sinO = 0.5 and C^^^ = 10 m^-. 
Increasing values of the ratio causes to increase \Acp{Afb)\ for Gj-^-^ > 0. With the increasing mass 
ratio of and ^4°, \Acp{Afb) \ can take large values. For Gj-^-^ < 0, Acp{Afb) is not sensitive to 
the mass ratio. 



9 



Now, we would like to summarize our results. 

• \Apb\ for the process under consideration is at the order of 10"^ and smaller compared to the 
SM one, for Cj^^ > 0. It can exceed the SM value (0.195) for Cj^^ < 0. Addition of NHB 
effects decreases its magnitude by 30% (slightly) for Cy-^-^ > (Cj-^'^ < 0). ApB is sensitive to 
the parameters sin6, .^^ and especially for C^-^^ > 0. Its magnitude decreases (increases) 
with increasing values of (^^)- 

• \Acp\ is at the order of 10~^. Addition of NHB effects decreases its magnitude by 50% 
(slightly) for C^'^^ > (C^^^^ < 0). It has the same sign in the restriction region C^^^ > 
and it can take both signs for C^^^ < 0. Acp is sensitive to the parameters sinO, espe- 
cially for Cy^-^ > 0. It decreases with increasing values of ^n,tt- The sensitivity of Acp to the 
ratio is weak. 

• Acp{Afb) is at the order of the magnitude of 1% for the intermediate values of sin9 for Cy-^-^ > 
0. It has the same sign in the restriction region Cj^^ > and it can take both signs for C^-^^ < 
0. Acp{Apb) is sensitive to the parameters and for Cj'^-^ > 0. It increases with 
increasing values of ^n,tt'> ^^en reach to 10%. Further the increasing values of the ratio 
causes to increase \Acp{Afb)\. 

Therefore, the experimental investigation of App and Acp and Acp{Afb) ensure a crucial test 
for new physics effects beyond the SM and also the sign of C^^-^. 
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A The operator basis 

The operator basis in the 2HDM (model III ) for our process is [ [T^ , ffT], H3] 

q=u,d,s,c,b 
q=u,d,s,c,b 

O5 = (sLa^tibla) (QRlBYQRf^), 
q=u,d,s,c,b 

q=u,d,s,c,b 

g 

Os = j^s^T^^a^,{m,R + msL)bf,g''>''' , 
g 

g 

= 4r^islb%)ifr) 

Qs = 4^2(^1^ r) E ^t^i) 

q=u,d,s,c,h 
q=u,d,s,c,b 
q=u,d,s,c,b 
q=u,d,s,c,b 

q=u,d,s,c,b 

QS = ^2(-'LCT^''b%) E ifR^^^Ql) 

q=u,d,s,c,b 

Q9 = ^2(''l^"'b'R) E ifL^,M 

q=u,d,s,c,b 

Qio = ]^2i'L^"'b''n) E i<lR^,>^Ql) 

q=u,d,s,c,b 
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where a and (5 are SIJ{?)) colour indices and T^^ and Q^^ are the field strength tensors of the elec- 
tromagnetic and strong interactions, respectively. Note that there are also flipped chirality partners of 
these operators, which can be obtained by interchanging L and R in the basis given above in model 
III. However, we do not present them here since corresponding Wilson coefficients are negligible. 

B The Initial values of the Wilson coefficients. 

The initial values of the Wilson coefficients for the relevant process in the SM are pT} ] 

Ct'imw) = 1, 
■^SM/ \ _ — 2xf —8xf — 5xf + 7xt 



4(xt - ly 2A{xt - 1; 



3 



^"^^^ - ~4(x,-l)^^^^*+ 8(x,-l)3 ' 

C-^(m^) = --^Bix,) + ' - ^^ Cix,) - D{x,) + I 
stn u\Y sm U\Y y 

C^o^'imw) = -4^iB{xt)-C{x,)) , 
sm a^r 

C^firnw) = z = l,..,10. (25) 
and for the additional part due to charged Higgs bosons are 

Ct.eimw) = 0, 

Cf(m^) = F.ivt) + XY F2iyt) , 

C^imw) = Y^ Civt) + XY G2iyt) , 

Ci^(miy) = Y'L^{yt) , (26) 



where 



Y - (^D ,?D ^ts\ 

^ — I ^N,bb "T ^N,sbT7' j ) 

nib \ ' VtbJ 

y = — Uht + iltc'^] , (27) 

mt \ ' ' Vt*J 



The NHB effects bring new operators and the corresponding Wilson coefficients read as 
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^Q2 ((^^,tt)^) 



^Q2 i^^fit) — 



Qi{zA,yt)05{yt) 



xyxt + 2za 2 Ij^^ ^^®5 ) 



-xyxtyti- 



dAn'^m^omt \Q5ixt) <di{zA,Xt) Qi{zA,Xt) 
{x - l)xt{yt/zA - 1) - (1 + x)yt) 



2pD tD 



xjyt + xtjl - yt/zA)) - 2yt 
QeQsizA) 



XfVt + 2y(a; - l)xtyt - za{xI + ©e) ^ ytjzj^ ^ ^ zaQ% 



(66 -{x- y){xt - yt))iQ3izA) + (x - y){xt - yt)zA) 

^QA^N,bb)- g4^2^2 



('AO 



r^H^atU n2n _ 9'^{iN,tt)'^'fnbmr ( Xt{l - 2y)yt , {-I + 2 cof? 6 w) {-I + x + y)yt 



cos^ OwQiivt) 



^ zh{Qi{zh, yt)xyt + cos^ Ow (-2x^(-l + a:t)|/|/| + xxtyyf - ©8^g)) \ 

cos2 6'w^ei (2:^7, 7/4)67 J' 



(28) 



-1 + 2cos^ gyt/) i/t _ xtyt _^ xtytjxy - zh) 
^biVt) QiizH^yt) 

Qb{.yt)zH 



cos2 Ow QiiVt) 

-1 + 2cos^ dw)ytZH 
cos^ ^14^67 



2xt In 



.0i(^i?,yt). 

4 _ g^mbin^xt / _ ^ , (-1 + 2x)xt , 2a;t(-l + (2 + xt)y) 



^1 ^ 1287r2m|,om2 1^ ^ 65(0;,) + y{l - Xt) 



+ 



65 (xt) 



4cos2 6li^(-l + X + y) + xt{ x + y) ^ Xt{x{xt {y - 2zh) - ^zh) + 2zh) 



cos^ 9wQi{xt) 



Ql{zH:Xt) 



yt((-l + x)xt2;ij + cos^ Owii^x - y)zH + Xt{2y{x - 1) - 2;if(2 - 3x - y)))) 



cos2 6'iv(63(^h) + - yt)zH) 



+2 {xt In 



0i(^ij,2;t) 



+ In 



xjyt - Xt)zH - Qsizn) 
(65 (xt) +y{l-Xt)ytZH 



32n^mlomtQi{zh, 1/4)65(2/*) 



(Q,{zh, yt){2y - 1) + 65(yt)(2x - l)zh) 



rho((cU ^2^ _ ^^,rr(C^,tt)^ {^b{yt)zh{yt - l)(x + y - 1) - 6i(^ft, yt){QM + ^t) iv. 



61(^^)65(^4) 



^fc€>5(yt) 
6i(zft) . 



^^^^,rr^^,ti"^6'^t /2(-l + (2 + Xt)?/) Xt{x - l){yt - Zh) 



65(Xt) 



+ 



x{xt{y - 2zh) - 4:Zh) + 2zh 

Ql{Zh,Xt) 



e'2{zh) 
(1 + x)ytZh 



+ 2 In 



^^65(3:4) 

_ei{zh,xt)_ 



xyxtyt + Zh{{x- y){xt - yt) - 65) 



^ 69 + ytZh{{x - y){xt - yt) - 65) (2a; - 1) ^ x{ytZh + xtjzh - yt)) - 2ytZh \ 
^ft 65(66 -{x- y){xt - yt)) 62(2;^) J ' 
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9^^N,rT^N,bb ( yXtVtixX^iyt - Zh) + Q(iZh{x - 2)) 



QAir'^mf^o \ 



+ 21n 



XtVt 



+ 2 In 



XtVtZh 



where 



ei(cj,A 

e5(A 

e, 



67 
eg 



= — (— 1 + y — ?/A)u; — a;(?/A + — cjA) 

= (a;* + 2/(1 - a;i))?/tt^ - xxtiyyt + (y* - l)u;) 

= Q2{'^,xt^ yt) 

= {xt{-l + y)- y)yti^ + xxtiyyt + ^^(-1 + yt)) 

= 1 — X + XLJ 

= X + A(l — x) 

= {xt + y{l - xt))yt + xxt{l - yt) 

= {yiVt - 1) - yt)zH + x{yyt + {yt - 1)zh) 

= yt{2x\l + xt){yt - 1) + xt{y{l - yt) + yt) + x{2{l -y + yt) 

+ xtil-2y{l-yt)-3yt))) 

= -x^{-l + X + y){-yt + x{2yt - l)){yt - Zh) - xtytZh{x{l + 2x) - 2y) 

+ VtiM^"^ - 2/(1 - a;)) + (1 + x)(x - 



(29) 



and 



Vt 



mH± 



rat 



mt 



Zh 



m 



Zh 



HO 



m 



za 



m 



AO 



The explicit forms of the functions Fi(2)(yf), Gi{2){yt), Hi{yt) and Li{yt) in eq.(|^) are given as 

yt{7-byt-^y^,) , y^ti^yt-2) 



Fiiyt) 
F2{yt) 
Gliyt) 
G2{yt) 
Hi{yt) 



72{yt - 1)3 + 12(y, - 1)4 
yt{5yt-3) , yt{-3yt + 2) ^^^^ 



+ 



12{yt - 1)2 6{yt - 1) 
yt{-yl + byt + 2) , 



-yt 



2^{yt - If 



Kyt - 1)' 



+ 



yt 



ytjyt - 3) 

%,-i)2 ' 2(y,-l)3 
1 — Asin'^Ow xyt 



siv?6 



w 



1 



- yt 



Liiyt) 



yt - 1 {yt- 1)2 
477/2 -7% + 38 3i/3-6|/, + 4 
18(2/t - 1)4 
1 



108(2/* - 1)=^ 



sin^Ow 8 



+ 



y* - 1 (i/i - 1)2 



(30) 
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Finally, the initial values of the coefficients in the model III are 



uir sm UW 

Cj™^(mH.) = , ^ = 5,...,10. (31) 

Here, we present Cq^ and Cq^ in terms of the Feynmann parameters x and y since the integrated 
results are extremely large. Using these initial values, we can calculate the coefficients Cf^^^^{n) 
and Cq^^^\ix) at any lower scale in the effective theory with five quarks, namely u, c, d, s, b similar 
to the SM case [0 ^ ^ 



The Wilson coefficients playing the essential role in this process are Cj^^^ (fi), Cl^^^ {fj), 
(j2HDM^^^^ Q'mDAi^^-^ (j-niDMi^^y p^j. completeness, in the following we give their explicit 

expressions. 

where the LO QCD corrected Wilson coefficient Cj^'^^^^^fi) is given by 

^LO,2HDM^^^ = 7^16/23^2™/ ^ (8/3) (r^^^/^S _ ^ 16/23 )^2i/DM(^^) 

+ C2'™^(mH/)EM'^^ , (32) 

i=l 

and T] = as{mw) /ctsifJ'), hi and are the numbers which appear during the evaluation [|T3|]. 

Cg-^-^ (fi) contains a perturbative part and a part coming from LD effects due to conversion of the 
real cc into lepton pair r+r^ : 

C^^if^) = Cr\f^) + Yresonis) , (33) 

where 
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ri2HDM 



+ h{z, s) (3Ci(/i) + C2{^^) + 3C3(/i) + C^{^^) + + Cel/i)) 

- 5) (4C3(/i) + 4C4(/X) + ?,CM + C6(/X)) 

- ^/i(0, s) {CM + 3Q(^)) + ^ (3C3(/x) + C4(/i) + 3C5(/i) + C6(/i)) 



(34) 



and 



Y 



a 



em Vi=ipi 



T ' T 



- my. + imyT 



(3Ci(/i) + C2(/x) + SCsifi) + C4(/i) + 3C,ifi) + C6(^)) . 
In eg .(1331), the functions h{u, s) axe given by 

h{u, s) 



m, 8 ^84 

— m m M H — X 

27 9 



9 /i 9 

-^(2 + a;)|l 
9 



(35) 



(36) 



In 



1-x-l 

2 arctan / , . 



— , for X = — < 1 



for a; = ^ > 1 



/i(0,s) 



8 8, mfe 4 4. 

m m s H — «7r 

27 9 /i 9 9 



(37) 



with u = The phenomenological parameter k in eq. ( |35| ) is taken as 2.3. In eqs. (30) and (|35|), 
the contributions of the coefficients Ci{fi), Cq{ii) are due to the operator mixing. 
Finally, the Wilson coefficients Cq^ (fi) and Cq^ (/i) are given by [|16p 



CQXf^) = V~''^''CQXmw),^ = 1,2. 

C The form factors for the decay B ^ K*l^l 

The structure functions appearing in eq. ( [TT| ) are 

-± 4C?"!^r., 

niB + rriK' q 



(38) 



A = 




Bi = 


r^ff 


B2 = 


r^ff 


Bs = 


r^ff 


C = 





A. 



mB + rriK 



ff^b 



T2 + 



-^3 , 



{A,-A,)+AC^,ff'^T, 



V 



niB + rriK* 
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n 


— 


D2 


= -c, 


n 


— n 


F 


= C'qi 


G 


= C'q2 



A, 



10 



1713+ TTLk* 



mi, 
2mK* 
rrib 



(39) 



We use the dependent expression which is calculated in the framework of light-cone QCD sum 
rules in p3| ] to calculate the hadronic form factors V, Ai, A2, Aq, Ti, T2 and T3: 

m 



1 - apK + bpiK)'' 



(40) 



where the values of parameters -F(O), ap and hp are listed in ^ 





m 


ap 


bp 


Ai 


0.34 ±0.05 


0.60 


-0.023 


A2 


0.28 ±0.04 


1.18 


0.281 


V 


0.46 ±0.07 


1.55 


0.575 


Ti 


0.19 ±0.03 


1.59 


0.615 


T2 


0.19 ±0.03 


0.49 


-0.241 


T3 


0.13 ±0.02 


1.20 


0.098 



Table 2: The values of parameters existing in eq. 
B K*. 



1) for the various form factors of the transition 
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Figure 1: ApB sl^sl function of sinO for m^o = 80 GeV without NHB effects. Here Aps is restricted 
in the region between solid (dashed) lines for C^'^'^ > (C^^^ < 0). Straight line corresponds to the 
SM contribution. 
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sind 



Figure 2: The same as Fig.|I], but for = 10 111^ and including NHB effects. 
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Figure 3: ApB as a function of ^^^^ for sin6 = 0.5 and 171^0 = 80 GeV. Here Aps is restricted in 
the region between solid (dashed) lines for Cy-^^ > (C^-^-^ < 0). Straight line corresponds to the SM 
contribution. 




Figure 4: ApB as a function of for sin9 — 0.5 and .C^,rr = 10 "^t- Here Afb is restricted in the 

region between solid (dashed) lines for C^-^^ > (Cy^^ < 0). Straight line corresponds to the SM 
contribution. 
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Figure 7: The same as Fig. || but for Acp. 
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Figure 8: The same as Fig. ^ but for Acp. 



24 




0.2 0.3 



0.4 0.5 0.6 0.7 
sm6 



0.8 



Figure 9: The same as Fig. |^but for Acp{Afb)- 
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Figure 10: The same as Fig. || but for Acp{Afb)- 
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Figure 11: The same as Fig. ^ but for Acp{Afb)- 
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